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Abstract 

We calculate the form factor f+{q'^) for i?-meson semileptonic decay in unquenched lattice QCD 
with 2+1 flavors of light sea quarks. We use Asqtad-improved staggered light quarks and a Fermilab 
bottom quark on gauge configurations generated by the MILC Collaboration. We simulate with 
several light quark masses and at two lattice spacings, and extrapolate to the physical quark mass 
and continuum limit using heavy-light meson staggered chiral perturbation theory. We then fit 
the lattice result for simultaneously with that measured by the BABAR experiment using 

a parameterization of the form factor shape in q'^ which relies only on analyticity and unitarity in 
order to determine the CKM matrix element |V^fe|. This approach reduces the total uncertainty 
in \ Vui,\ by combining the lattice and experimental information in an optimal, model-independent 
manner. We find a value of |Kb| x 10^ = 3.38 ± 0.36. 

PACS numbers: 12.15.Hh, 12.38. Gc, 13.20.He 
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I. INTRODUCTION 



The semileptonic decay B iriu is a sensitive probe of the heavy-to-hght quark-flavor 
changing h ^ u transition. When combined with an experimental measurement of the 
differential decay rate, a precise QCD determination of the B — > ttIv form factor allows a 
clean determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vub\ with 
all sources of systematic uncertainty under control. In the Standard Model, the differential 
decay rate for this process is 

where q = Pb — Pn the momentum transferred from the i?- meson to the outgoing leptons. 
The form factor, parameterizes the hadronic contribution to the weak decay, and 

must be calculated nonperturbatively from first principles using lattice QCD. 

A precise knowledge of CKM matrix elements such as \Vub\ is important not only be- 
cause they are fundamental parameters of the Standard Model, but because inconsistencies 
between independent determinations of the CKM matrix elements and CP-violating phase 
would provide evidence for new physics. Although the Standard Model has been amazingly 
successful in describing the outcome of most particle physics experiments to date, it cannot 
account for gravity, dark matter and dark energy, or the large matter-antimatter asymmetry 
of the universe. Thus we know that it is incomplete, and expect new physics to affect the 
quark-flavor sector to some degree, although we do not know a priori what experimental 
and theoretical precision will be needed to observe it. 

The determination of \Vub\ from B — ^ iriu semileptonic decay relies upon the assump- 
tion that, because the leading Standard Model decay amplitude is mediated by tree-level 
VT-boson exchange, it will not be signiflcantly affected by new physics at the current level 
of achievable precision. Recently, however, hints of new physics have appeared in various 
regions of the heavy-quark flavor sector such as CP-asymmetries in P ^ Kir [1], constraints 
on sin(2/3) from AP = 2 neutral meson mixing and 1-loop penguin-induced decays [2], and 
the phase of the P^-mixing amplitude [3-5]. The unexpected inconsistency most relevant 
to our new lattice QCD calculation of the P — niu form factor and \Vub\ is the current 
"/ds puzzle" [6]. The HPQCD Collaboration's lattice QCD calculation of the P^-meson 
leptonic decay constant /^^ [7] disagrees with the latest results from the Belle, BABAR, 
and CLEO experiments [8-12] at the 3-cr level, although HPQCD's determinations of the 



masses m£)+ and ttids and the decay constants /tt, fx, and //)+ all agree quite well with ex- 
perimental measurements [13, 14]. Furthermore, because the significance of the discrepancy 
is dominated by the statistical experimental uncertainties, it cannot easily be explained by 
an underestimate of the theoretical uncertainties. Additional lattice QCD calculations of 
Jd^ are needed to either confirm or reduce the inconsistency. If the disagreement holds up, 
however, it is evidence for a large new physics contribution to a tree-level Standard Model 
process at the few percent-level. Therefore, although B rriu semileptonic decay provides 
a theoretically clean method for determining \Vub\ within the framework of the Standard 
Model, we should keep in mind that new physics could appear in 6 — > m transitions. 

Understanding and controlling all sources of systematic uncertainty in lattice QCD cal- 
culations of hadronic weak matrix elements is essential in order to allow accurate determi- 
nations of Standard Model parameters and reliable searches for new physics. The hadronic 
amplitudes for B Tiiu, in particular, can be calculated accurately using current lattice 
QCD methods because the decay process is "gold plated", i.e., there is only a single stable 
hadron in both the initial and final states. Lattice calculations with staggered quarks al- 
low for realistic QCD simulations with dynamical quarks as light as m^/lO, multiple lattice 
spacings, large physical volumes, and high statistics. The resulting simulations of many 
light-light and heavy-light meson quantities with dynamical staggered quarks are in excel- 
lent numerical agreement with experimental results [15]. This includes both postdictions, 
such as the pion decay constant [16], and predictions, as in the case of the B^ meson mass 
[17]. Such successes show that the systematic uncertainties in these lattice QCD calculations 
are under control, and give confidence that additional calculations using the same methods 
are reliable. 

The publicly available MILC gauge configurations with three fiavors of improved stag- 
gered quarks [18] that have enabled these precise lattice calculations make use of the 
"fourth-root" procedure for removing the undesired four-fold degeneracy of staggered lat- 
tice fermions. Although this procedure has not been rigorously proven correct, Shamir uses 
plausible assumptions to argue that the continuum limit of the rooted theory is in the same 
universality class as QCD [19, 20]. The rooting procedure leads to violations of unitar- 
ity that vanish in the continuum limit; both theoretical arguments [21, 22] and numerical 
simulations [23-25], however, show that the unitarity- violating lattice artifacts in the pseudo- 
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Goldstone boson sector can be described and hence removed using rooted Staggered Chiral 
Perturbation Theory (rS^PT), the low-energy effective description of the rooted staggered 
lattice theory [26-28]. Given the wealth of numerical and analytical evidence supporting 
the validity of the rooting procedure, most of which is reviewed in Refs. [29-31], we work 
under the plausible assumption that the continuum limit of the rooted staggered theory is 
QCD. We note, however, that it is important to have crosschecks of lattice calculations of 
phenomenologically-important quantities using a variety of fermion formulations, since they 
all have different sources of systematic uncertainty. 

Both existing unquenched lattice calculations of the B — > iriu form factor use the MILC 
configurations. When combined with the Heavy Flavor Averaging Group's latest determi- 
nation of the experimental decay rate from ICHEP 2008 [32], they yield the following values 
for \Vub\: 

\Vub\ X 10^ = 3.40 ± 0.20t[J;i HPQCD [33], (2) 

\Vub\ X 10^ = 3.62 ± 0.221°;^^ Fermilab-MILC [34], (3) 

where the errors are experimental and theoretical, respectively. Both analyses primarily rely 
upon data generated at a "coarse" lattice spacing of a ~ 0.12 fm, and use a smaller amount 
of "fine" data at a ~ 0.09 fm to check the estimate of discretization errors. Neither is 
therefore able to extrapolate the B —* niu form factor to the continuum {a 0). The most 
significant difference in the two calculations is their use of different lattice formulations for 
the bottom quarks. The HPQCD Collaboration [33] uses nonrelativistic (NRQCD) heavy 
quarks [35], whereas we use relativistic clover quarks with the Fermilab interpretation [36] 
via heavy quark effective theory (HQET) [37-39]. Both methods work quite well for heavy 
bottom quarks. The Fermilab treatment, however, has the advantage that it can also be 
applied to charm quarks; we can therefore use the same method for other semileptonic 
form factors such as -D ^ niu, D Kiu, and B D*iv [40, 41]. The two unquenched 
lattice calculations of the B iriu form factor, which have largely independent sources 
of systematic uncertainty, nevertheless lead to consistent values of \Vub\ with similar total 
errors of ~ 15%. 

In this paper we present a new model-independent unquenched lattice QCD calculation 
of the B -iriu semileptonic form factor and \Vub\- Our work builds upon the previous 
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Fermilab-MILC calculation and improves upon it in several ways. We now include data on 
both the coarse and fine MILC lattices, and can therefore take the a —>■ limit of our data 
which is generated at nonzero lattice spacing. We also have additional statistics on a subset 
of the coarse ensembles. The most important improvements, however, are in the analysis 
procedures. 

We have removed all model-dependent assumptions about the shape in of the form 
factor from the current analysis. Our result is therefore theoretically cleaner and more 
reliable than those of previous lattice QCD calculations. The first refinement over previous 
unquenched lattice B txIv form factor calculations is in the treatment of the chiral and 
continuum extrapolations. We simultaneously extrapolate to physical quark masses and zero 
lattice spacing and interpolate in the momentum transfer by performing a single fit to our 
entire data set (all values of m^, a, and q^) guided by functional forms derived in heavy-light 
meson staggered chiral perturbation theory (HMSxPT) [42] . We thereby extract the physical 
form factor f+{q^) in a controlled manner without introducing a particular ansatz for the 
form factor's dependence. The second refinement over previous unquenched B niu 
lattice form factor calculations is in the combination of the lattice form factor result and 
experimental data for the decay rate to determine the CKM matrix element \Vub\- We fit our 
lattice numerical Monte Carlo data and the 12-bin BABAR experimental data [43] together 
to the model- independent "2;-expansion" of the form factor given in Ref. [44], in which the 
form factor is described by a power series in a small quantity z with the sum of the squares of 
the series coefficients bounded by unitarity constraints. We leave the relative normalization 
factor, \Vub\, as a free parameter to be determined by the fit, thereby extracting \Vub\ in 
an optimal, model-independent way. Others have also fit lattice and experimental results 
together using different, equally- valid, parameterizations [45, 46]. This work, however, is the 
first to use the full correlation matrices, derived directly from the data, for both the lattice 
calculation and experimental measurement. 

This paper is organized as follows. In Sec. II we describe the details of our numerical 
simulations. We discuss the gluon, light-quark, and heavy-quark lattice actions, and present 
the parameters used, such as the quark masses and lattice spacings. We then define the 
matrix elements needed to calculate the semileptonic form factors and discuss the method for 
matching the lattice heavy-light current to the continuum. Next we describe our analysis for 
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determining the form factors in Sec. III. This is a three-step procedure. We first fit pion and 
5-meson 2-point correlation functions to extract the meson masses. We then fit the B ^ ir 
3-point function, using the masses and amphtudes from the 2-point fits as input, to extract 
the lattice form factors at each value of the light quark mass and lattice spacing. Finally, 
we extrapolate the results at unphysical quark masses and nonzero lattice spacing to the 
physical light quark masses and zero lattice spacing using HMSxPT. In Sec. IV we estimate 
the contributions of the various systematic uncertainties to the form factors, discussing each 
item in the error budget separately. We then present the final result for f+{q^) with a 
detailed breakdown of the error by source in each bin. We combine our result for the 
form factor with experimental data from the BABAR Collaboration to determine the CKM 
matrix element \Vub\ in Sec. V. We also define the model-independent description of the 
form factor shape that we use in the fit and discuss alternative parameterizations of the form 
factor. Finally, in Sec. VI we compare our results with those of previous unquenched lattice 
calculations. We also compare our determination of \Vub\ with inclusive determinations and 
to the preferred values from the global CKM unitarity triangle analysis. We conclude by 
discussing the prospects for improvements in our calculation and its impact on searches for 
new physics in the quark flavor sector. 

II. LATTICE CALCULATION 

In this section we describe the details of our numerical lattice simulations. We first present 
the actions and parameters used for the light (up, down, strange) and heavy (bottom) quarks 
in Sec. II A. We then define the procedure for constructing lattice correlation functions 
with both staggered light and Wilson heavy quarks in Sec. II B. Finally, in Sec. II C, we 
show how to match the lattice heavy-light vector currents to the continuum with a mostly 
nonperturbative method, so that lattice perturbation theory is only needed to estimate a 
small correction. 

A. Actions and Parameters 

We use the ensembles of lattice gauge fields generated by the MILC Collaboration and 
described in Ref. [18] at two lattice spacings, a ^ 0.12 and 0.09 fm, in our numerical lattice 
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TABLE I: Lattice simulation parameters. The cokimns from left to right are the approximate 
lattice spacing in fm, the bare light quark masses ami/anis, the linear spatial dimension of the 
lattice in fm, the dimensionless factor m^^L (corresponding to the taste-pseudoscalar pion composed 
of light sea quarks), the dimensions of the lattice in lattice units, the number of configurations used 
for this analysis, the clover term csw and bare k value used to generate the bottom quark, and 
the improvement coefficient used to rotate the bottom quark field in the 6 — > u vector current. 



a(fm) 


ami 1 iiT^s 


L(fm) 




Volume 


# Configs. 


csw Kb 


di 


0.09 


0.0062/0.031 


2.4 


4.1 


28^ X 96 


557 


1.476 0.0923 


0.09474 


0.09 


0.0124/0.031 


2.4 


5.8 


28^ X 96 


518 


1.476 0.0923 


0.09469 


0.12 


0.005/0.05 


2.9 


3.8 


24^ X 64 


529 


1.72 0.086 


0.09372 


0.12 


0.007/0.05 


2.4 


3.8 


20^ X 64 


836 


1.72 0.086 


0.09372 


0.12 


0.01/0.05 


2.4 


4.5 


20^ X 64 


592 


1.72 0.086 


0.09384 


0.12 


0.02/0.05 


2.4 


6.2 


20^ X 64 


460 


1.72 0.086 


0.09368 



simulations of the B Triu semileptonic form factor. The ensembles include the effects 
of three dynamical staggered quarks — two degenerate light quarks with masses ranging 
from ms/8 - ms/2 and one heavier quark tuned to within 10-30% of the physical strange 
quark mass. This allows us to perform a controlled extrapolation to both the continuum 
and the physical average u-d quark mass. The physical lattice volumes are all sufficiently 
large [rrij^L > 4) to ensure that effects due to the finite spatial extent remain small. 

For each independent ensemble we compute the light valence quark in the 2-point and 
3-point correlation functions only at the same mass, mi, as the light quark in the sea sector. 
Thus all of our simulations are at the "full QCD" point. Note, however, that we still have 
many correlated data points on each ensemble because of the multiple pion energies. Table I 
shows the combinations of lattice spacings, lattice volumes, and quark masses used in our 
calculation. 

For bottom quarks in 2-point and 3-point correlation functions we use the Sheikholeslami- 
Wohlert (SW) "clover" action [47] with the Fermilab interpretation via HQET [36, 37], which 
is well-suited for heavy quarks, even when amq > 1. Because the spin- flavor symmetry 
of heavy quark systems is respected by the lattice regulator, the expansion in l/rriQ of 
the heavy-quark lattice action has the same form as the l/mq expansion of the heavy- 
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quark part of the QCD action. Discretization effects in the lattice heavy-quark action 
are therefore parameterized order-by-order in the heavy-quark expansion by deviations of 
effective operator coefficients from their values in continuum QCD. Thus, in principle, the 
lattice heavy-quark action can be improved to arbitrarily high orders in HQET by tuning a 
sufficiently large number of parameters in the lattice action. In practice, we tune the hopping 
parameter, k, and the clover coefficient, csw, of the SW action, to remove discretization 
effects through next-to-leading order, Oll/mq), in the heavy-quark expansion. 

The SW action includes a dimension-five interaction with a coupling csw fhat must be 
adjusted to normalize the heavy quark's chromomagnetic moment correctly [36]. In our 
calculation we set the value of csw = Uq^, as suggested by tadpole-improved, tree- level 
perturbation theory [48]. We determine the value of Uq either from the plaquette (a ~ 0.09 
fm) or from the Landau link (a ^ 0.12 fm). The tadpole-improved bare quark mass for SW 
quarks is given by 



such that tuning the parameter k to the critical quark hopping parameter K^it leads to a 
massless pion. Before generating the correlation functions needed for the B iriu form 
factor, we compute the spin-averaged Bg kinetic mass on a subset of the available ensembles 
in order to tune the bare k value for bottom (and hence the corresponding bare quark mass) 
to its physical value [36]. We then use the tuned value of Hh for the B niu form- factor 
production runs. Table I shows the values of the clover coefficient and tuned Kb used in our 
calculation. 

In order to take advantage of the improved action in the calculation of the B niu 
form factor, we must also improve the flavor-changing vector current to the same order in 
the heavy-quark expansion. We remove errors of 0{l/mQ) in the vector current by rotating 
the heavy-quark field used in the matrix element calculation as 



where Diat is the symmetric, nearest-neighbor, covariant difference operator. We set di to 
its tadpole- improved tree-level value of [36] 




(4) 




(5) 




(6) 
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The values of the rotation parameter used in our calculation are given in Table I. 

In order to convert dimensionful quantities determined in our lattice simulations into 
physical units, we need to know the value of the lattice spacing, a, which we find by com- 
puting a physical quantity that can be compared directly with experiment. We first deter- 
mine the relative lattice scale by calculating the ratio ri/a on each ensemble, where ri is 
related to the force between static quarks, r\F{ri) = 1.0 [49]. These ri/a estimates are then 
smoothed by fitting to a smooth function of the gauge coupling and quark masses. This 
scale-setting method has the advantage that the ratio ri/a can be determined precisely from 
a fit to the static quark potential [50, 51]. We convert all of our data from lattice spacing 
units into ri units before performing any chiral fits in order to account for slight differences 
in the value of the lattice spacing between ensembles. In this work we use the value of 
^phys _ 0.3108(15) (1; 79) obtained by combining a recent lattice determination of ri/^ [52] 
with the PDG value of = 130.7 ± 0.1 ± 0.36 MeV [53] to convert lattice results from ri 
units to physical units. 



B. Heavy-light meson correlation functions 

The B irlu semileptonic form factors parameterize the hadronic matrix element of the 
b —>■ u quark flavor-changing vector current = iu'-f'^b: 

where is the momentum transferred to the outgoing lepton pair. For calculations on the 
lattice and in HQET, it is more convenient to write the matrix element as [54] 

where = p'^/rriB is the four- velocity of the 5- meson, = — {pt^-v)v^ is the component 
of the pion momentum orthogonal to f, and E^^ = p.,^ ■ v = (m^ + — q^)/{2mB) is the 
energy of the pion in the 5-meson rest frame {pb = 0). In this frame the form factors f\\ {E^^) 
and /±(-E'7r) are directly proportional to the hadronic matrix elements: 
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We therefore first calculate the hadronic matrix elements in Eqs. (9) and (10) in the rest 
frame of the 5-meson to obtain f\\{ETr) and f±{ETr), and then extract the standard form 
factors foiq"^) and f+{q^) using the following relations: 

/o(g') = [{ms-E^)f\\{E^) + {El-ml)f^{E^)\, (11) 

U{q') = [/||(E.) + {ms - E^)h{E^)] . (12) 

These relations automatically satisfy the kinematic constraint /+(0) = /o(0). 

The 2-point and 3-point correlation functions needed to extract the lattice matrix element 
for B iriu decay are 

CUt;P.) = X]e*^^-^^(a(0,0)Ot(t,f)), (13) 

X 

c^it) = J2(^Bio,6)oUt,x)), (14) 

X 

Cl-;^{t,T-p^) = Y^^e^^^^-y{O^{0,6)V,{t,y)Ol{T,x)), (15) 

x,y 

where Ob and (9^ are interpolating operators for the 5- meson and pion, respectively, and 
Vn is the heavy-light vector current on the lattice. 

In practice, to construct the heavy-light bilinears we must combine a staggered light 
quark, which is a 1-component spinor, with a 4-component Wilson-type bottom quark; we 
do so using the method established by Wingate et al. in Ref. [55]. For the B meson we use 
a mixed-action interpolating operator: 

OBd^) = i^^{xhlpn^E{x)x{x), (16) 

where a, j3 are spin indices and f2(x) = To^Ti ^72^73^- The fields ip and x a^'s the 4-component 
clover quark field and 1-component staggered field, respectively. Based on the transformation 
properties of Ob,-e{x) under shifts by one lattice spacing, S plays the role of a (fermionic) 
taste index [31, 55]. Once Ob,-e{x) is summed over 2^ hypercubes in the correlation functions 
that we compute, S also takes on the role of a taste degree of freedom, in the sense of 
Refs. [56, 57]. Because the heavy quark field ^Ja^x) is slowly varying over a hypercube, it 
does not affect the construction of Refs. [56, 57]. 

For the pion we use the local pseudoscalar interpolating operator, 

0^{x)=e{x)x{x)x{x), (17) 
11 



where e{x) = (— 1) (^1+^2 +X3 +0:4) _ take the vector current to be 

Vi{x) = ^,(x)7^^fi/5E(x)x(a:), (18) 

where \l/ is the rotated heavy-quark field given in Eq. (5). When forming and 
C^^'^(t, T; Ptt), we sum over the taste index. This yields the same correlation functions, 
with respect to taste, as in Ref. [55]. Our principal difference with Ref. [55] is to use 4- 
component heavy quarks instead of 2-component non- relativist ic quarks, and to derive the 
correlators in the staggered formalism, without the introduction of naive fermions. 

We work in the rest frame of the 5-meson, so only the pions carry momentum. We 
compute both the 2-point function C^it^pT^) and the 3-point function C^^'^{t,T;pTr) at 
discrete values of the momenta p; = 27r(0, 0, 0)/L, 27r(l, 0, 0)/L, 27r(l, 1, 0)/L, 27r(l, 1, 1)/L, 
and 27r(2,0,0)/L allowed by the finite spatial lattice volume. We use only data through 
momentum j4 = 27r(l, 1, however, because the statistical errors in the correlators 
increase significantly with momentum. 

We use a local source for the pions throughout the calculation, while we smear 
the 5-meson wavefunction in both the 2-point function C^(t) and the 3-point function 

dndt^x) = J2Siy)'^ait,x + yhlpnp^it,x)xit,x), (19) 
y 

where S{y) is the spatial smearing function. This reduces contamination from heavier excited 
states and allows a better determination of the desired ground state amplitude. In our study 
of choices for how to smear the S-meson, we found that a wall source, S{y) = 1, worked 
extremely well for suppressing excited state contamination, but at the cost of large statistical 
errors in the 2-point and 3-point correlation functions. In contrast, use of a IS wavefunction, 
S{y) = exp(— /i|y|), optimized to have good overlap with the charmonium ground state led to 
smaller statistical errors at the cost of undesirably large excited state contributions to the 3- 
point function that would make it difficult to extract the ground state amplitude. In order to 
achieve a balance between small statistical errors and minimal excited state contamination, 
we tune the coefficient of the exponential in the IS wavefunction to the smallest value {i.e., 
the widest smearing) for which the 5-meson 2-point effective mass is still well-behaved; we 
find a value of afi = 0.20 for the coarse ensembles. We note that our determination of the 
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optimal S-meson smearing function is consistent with the theoretical expectation that the 
i?-meson wavefunction should be wider than the charmonium wavefunction. 

For the calculation of the 3-point function, we fix the location of the pion source at tj = 
and the location of the i?-meson sink at tj = T, and vary the position of the operator over 
all times t in between. If the source-sink separation is too small then the entire time range 
< t < T is contaminated by excited states, but if the source-sink separation is too large 
then the correlation function becomes extremely noisy. In practice, we set the sink time to 
T = 16 on the coarse lattices; we have checked, however, that our results using this choice 
are consistent with those determined from using T = 12 and T = 20. On the fine lattices we 
scale the source sink separation by the approximate ratio of the lattice spacings, ctfine/'^coarse, 
and use T = 24. 

In order to minimize the statistical errors given the available number of configurations 
in each ensemble, we compute the necessary 2-point and 3-point correlations not only with 
a source time of tj = 0, but also with source times of ti = nt/4:,nt/2, and Srit/A {nt is the 
temporal extent of the lattice) and the sink time T shifted accordingly. We then average 
the results from the four source times; this effectively increases our statistics by a factor of 
four. 



C. Heavy-light current renormalization 

In order to recover the desired continuum matrix element, the lattice amplitude must be 
multiplied by the appropriate renormalization factor Zy^: 

{7r\V,\B) = Zlx{n\V,\B), (20) 

where and are the lattice and continuum b ^ u vector currents, respectively. This 
removes the dominant discretization errors from the lattice current operator. In terms of 
the form factors, Eq. (20) can be rewritten as 

/|l = Z^f X (21) 
h = Z'^^ X ff, (22) 

where explicit expressions relating /||^* and to correlation functions are given in 
Eqs. (40) and (41). 

13 



In this work, we calculate Zy^ via the mostly nonperturbative method used in the earlier 
quenched Fermilab calculation [54]. We first rewrite Zy^ as 

Z'i=p'^,\fz^. (23) 

The flavor- conserving renormalization factors Zy and Zy account for most of the value of 
Zy [38]. They can be determined from standard heavy-light meson charge normalization 
conditions: 

Z^^ X = 1, (24) 

Z{* X {B\V'''''\B) = 1, (25) 

where the light-light and heavy-heavy lattice vector currents are given by 

V^iix) = xKxMx)i^Y^^n{x)p^,x{x), (26) 
V'^'^ix) = %M^'^^^,p{x), (27) 

respectively. In order to reduce the statistical errors in Zy, we compute the lattice matrix 
element {D\V'''''^\D) using a clover charm quark as the spectator in the 3-point correlation 
function. We eliminate contamination from staggered oscillating states in the determination 
of Zy by using a clover strange quark for the spectator in the 3-point correlation function 
{B\V'^^'^\B). Once Zy and Zy have been determined nonperturbatively, the remaining 
correction factor in Eq. (23), py^, is expected to be close to unity because most of the 
radiative corrections, including contributions from tadpole graphs, cancel in the ratio [38]. 
We therefore estimate py^ from 1-loop lattice perturbation theory [48]. 

The matching factor py^ has been calculated by a subset of the present authors, and a 
separate publication describing the details is in preparation [58]. The corrections to py^ can 
be expressed as a perturbative series expansion in powers of the strong coupling: 

p'^^ = l + A7rav{q*)Pv^'^ + 0{a'y), (28) 

where av{q*) is the renormalized coupling constant in the l^-scheme and is determined from 
the static quark potential with the same procedure as is used in Ref. [59]. The scale q*, 
which should be the size of a typical gluon loop momentum, is computed via an extension of 
the methods outlined by Brodsky, Lepage, and Mackenzie [48, 60] and Hornbostel, Lepage, 
and Morningstar [61]. The value of q* ranges from 2.0-4.5 GeV for the parameters used 

14 



in our simulations. The 1-loop coefficient, py^ , and higher moments are calculated using 
automated perturbation theory and numerical integration as described in Refs. [62, 63]. We 
find that the perturbative corrections to matrix elements of the temporal vector current, 
Vo, are less than a percent, while the corrections to matrix elements of the spatial vector 
current, Vi, are 3-4%. 

III. ANALYSIS 

In this section, we describe the three-step analysis procedure used to calculate the B — »• 
Tiiv semileptonic form factor, In the first subsection. Sec. IIIA, we describe how we 

fit the pion and i?-meson 2-point correlation functions in order to determine the pion energies 
and 5-meson mass. We use both of these quantities in the later determination of the lattice 
form factors f\\{ET^) and f±{ET^). Next, in Sec. IIIB, we construct a useful ratio of the 3- 
point correlation function (7r|y|i?) to the 2-point functions. We then fit this ratio to a simple 
plateau ansatz to extract the desired form factors. Finally, in Sec. Ill C, we extrapolate the 
form factors calculated at unphysically heavy quark masses and finite lattice spacing to 
the physical light quark masses and zero lattice spacing using next-to-leading order (NLO) 
HMSxPT expressions extended with next-to-next-to-leading order (NNLO) analytic terms. 
(We perform a simultaneous extrapolation in rriq and a and interpolation in Et^.) We then 
take the appropriate linear combination of fwiEj^) and f±{ET^) to determine the desired form 
factor, /+(g^), with statistical errors. 

A. Two-point correlator fits 

The pion and S-meson 2-point correlators obey the following functional forms: 



In the above expressions, terms with odd m contain the prefactor (—1)*. This leads to visible 
oscillations in time in the meson propagators; such behavior arises with staggered quarks 
because the parity operator is a composition of spatial inversion and translation through 




(29) 



(30) 
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one timeslice [64, 65]. The contributions of the opposite-parity oscillating states are found 
to be significant throughout the entire time range and must therefore be included in fits to 
extract the desired ground state energy. 

Because the statistical errors in the pion energies and 5-meson mass contribute very little 
to the total statistical error in the B Triu form factor, we use a simple procedure to fit 
the 2-point functions. Although this does not optimize the determinations of E^^ and rriB, it 
is sufficient for the purpose of this analysis. We first select a fit range, tmin^^max, that allows 
a good correlated, unconstrained fit including only contributions from the ground state and 
its opposite-parity partner. We then reduce tmin by one timeslice and redo the fit. If the 
correlated confidence level is too low (< 10%), we increase the number of states and try the 
fit again with the same time range. Otherwise, if the fit is good, we reduce t^[^ by one more 
timeslice and repeat the fit. We repeat this procedure until we can no longer get a good fit 
without using a large number (greater than 4) of states. We note that, by including only as 
many states as the data can determine, we minimize the possibility of spurious solutions in 
which the fitter exchanges the ground state with one of the same-parity excited states. We 
have, however, checked that this method yields the same results within statistical errors as 
a constrained fit that includes up to three or four pairs of states. 

Figure 1 shows examples of both vs. tmin (left plot) and ttib vs. tmin (right plot) 
on the ami/anis = 0.02/0.05 coarse ensemble, which has the largest light quark mass of 
the coarse ensembles. The masses are stable as tmin is reduced, and the statistical errors in 
rriB become smaller as additional timeslices are added to the fit. The statistical errors are 
determined by performing a separate fit to 500 bootstrap ensembles; each fit uses the full 
single elimination jackknife correlation matrix which is remade before every fit. The size of 
the statistical errors does not change when the number of bootstrap ensembles is increased 
by factors of two or four. We select the time range to use in the B iriu analysis based 
on several criteria: a good correlated confidence level, relatively symmetric upper and lower 
bootstrap error bars, no 5-a or greater outliers in the bootstrap distribution, and no sign 
of excited state contamination. The red (filled) data points in Fig. 1 mark the chosen fit 
ranges for the ensemble in the example plots. Figures 2 and 3 show the corresponding pion 
and 5-meson correlator fits, respectively, which go through the data points (shown with 
jackknife errors) quite well. 

The gauge configurations have been recorded every six trajectories, and the remaining 
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a = 0.12 fm; arni/arris = 0.02/0.05; tmax = 20 



a = 0.12 fm; ami/ arris = 0.02/0.05; fmax = 15 
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FIG. 1: Pion mass (left plot) and i?-meson mass (right plot) versus minimum timeslice in 2-point 
correlator fit. The red (filled) data points show the fit ranges selected for use in the B — > irlv form 
factor analysis. 
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FIG. 2: Pion correlator fit corresponding to the red data point in the left-hand graph of Fig. 1. 
The left plot shows the fit (red line) to the zero-momentum pion propagator on a log scale, while 
the right plot shows the deviation of the fit from the data point for each timeslice. On both plots 
the dashed vertical line indicates tmin- Single elimination jackknife statistical errors are shown. 



autocorrelations between consecutive configurations cannot be neglected. We address this 
by averaging a block of successive configurations together before calculating the correlation 
matrix and performing the fit. We determine the optimal block size by increasing the 
number of configurations in a block until the single elimination jackknife statistical error 
in the correlator data remains constant within errors. This is shown for a representative 
timeslice of the pion propagator on a coarse ensemble in Fig. 4. We find that it is necessary 
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a = 0.12 fm; ami/ am, = 0.02/0.05 



a = 0.12 fm; ami/am, = 0.02/0.05 




FIG. 3: S-meson correlator fit corresponding to tlie red data point in tlie right-liand grapli of 
Fig. 1. The left plot shows the fit (red line) to the i?-meson propagator on a log scale, while the 
right plot shows the deviation of the fit from the data point for each timeslice. On both plots the 
dashed vertical line indicates tmin- Single elimination jackknife statistical errors are shown. 
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FIG. 4: Single-elimination jackknife error versus block size in the zero-momentum pion propagator 
at t = 6. The statistical errors in the errors are calculated with an additional single elimination 
jackknife loop. The red line is an average of the errors for block sizes 5-8 and is only to make it 
easier to see that the statistical error plateaus after a block size of 5; it is not used in the form 
factor analysis. 



to use a block size of 5 on the coarse ensembles and 8 on the fine ensembles, and we use 
tliese values for the rest of the form factor analysis. We note that the size of the statistical 
errors that arises from blocking by 5 on the coarse ensemble is consistent with that estimated 
based on a calculation of the integrated autocorrelation time. 
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a = 0.12 fm; ami/ am, = 0.02/0.05 



a = 0.12 fm; ami /am, = 0.02/0.05 
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FIG. 5: Comparison of pion energy E^^ (left plot) and amplitude Z-,^ (right plot) with the prediction 
of the continuum dispersion relation. We also show a power-counting estimate for the size of 
momentum-dependent discretization errors, which are of 0{as{a\'PT^\)'^), as dashed black lines. 

The pion energy E^^ that is extracted from fitting the 2-point function, (t;p7r), should 
satisfy the dispersion relation E^. = |p7rP+^^ in the continuum limit due to the restoration of 
rotational symmetry. Similarly, the pion amplitude, = |(0|(9^|7r)|, should be independent 
of as a — s> 0. As shown in Fig. 5, our results are consistent with these continuum relations 
within statistical errors.^ We therefore replace the pion energy E^^ by \/j^5iT^~-fmJ when 
calculating the lattice form factors f\\{E.,^) and fxiE^^) in order to reduce the total statistical 
uncertainty. The pion amplitude drops out of the form factor calculation, however, because 
we take suitable ratios of 3-point to 2-point correlators. 



B. Three-point correlator fits 

The B ^ 71 3-point correlator obeys the following functional form: 



where 



''B 



^ As this analysis was being completed we generated data with four times the statistics on the ami Jams = 
0.02/0.05 coarse ensemble. In order to make the comparison to the continuum expectation clearer, we use 
the higher statistics data in Fig. 5. 
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Writing out the first four terms of C^^'^{t,T) makes the behavior of the 3-point correlator 
function of both t and T more transparent: 



As in the case of the pion and 5-meson propagators, the leading contributions from the 
opposite-parity excited states (the Ajji^ and A^^ terms) change sign when t ^ t + 1; these 
produce visible oscillations in the correlation function along the time direction. The sub- 
leading contribution from the opposite-parity excited states (the Aj^ term), however, only 
changes sign when the source-sink separation is varied, e.g., T — T+ 1; this contribution is 
not as clearly visible in the data as those that oscillate with the time slice t. 

The lattice form factors are related to the ground-state amplitude of the 3-point function 
C3^^-(t,T) as follows: 

fir - ^s" (^) (34) 

;...^.««(?^)_L, ,35) 

where, as before, Z.,, = |(0|(97r|7r)| and Zb = |(0|Cb|-B)|. The pion and 5-meson energies 
and amplitudes are known from the 2-point fits described in the previous subsection. Thus, 
the goal is to determine the 3-point amplitude A^^ for /i along both the spatial and temporal 
directions. 

In principle, the easiest way to determine the coefficient A^^ is to divide the 3-point 
function C^^'^{t,T) by the appropriate 2-point functions and fit to a constant (plateau) 
ansatz in a region of time slices ^ t ^ T that are sufficiently far from both the pion and B- 
meson sources, such that excited state contamination can be neglected. In practice, however, 
oscillating excited-state contributions are significant throughout the interval between the 
pion and 5-meson, so our raw correlator data cannot be fit to such a simple function. 
Therefore we construct an average correlator in which the oscillations are reduced before 
performing any fits. This method for determining the form factors requires knowledge of E^^ 
and rriB', we use the values determined in the 2-point fits described in the previous subsection 
and propagate the bootstrap uncertainties in order to properly account for correlations. 

The final ratio of correlators used to determine A^^ entails several pieces. To begin 
consider the carefully constructed average of the value of the 5-meson propagator at time 
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slice t with that at t + 1: 



(0), 
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2m 



(0)' 
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cgt) ^ + 

t g-m^"^(t+l) 
7/2 



(0) . 
-m.g't 



2m 



(1) 



(36) 



where Am^ = m^"* —m^^. By removing the leading exponential behavior from the correlator 
before taking the average we suppress the leading oscillating contribution by a factor of the 
mass-splitting AmB/2 while leaving the desired ground state amplitude unaffected. Note 
also that, while this procedure affects the size of the excited state amplitudes, it does not 
alter the functional form of the correlator, nor does it alter the energies in the exponentials. 
Therefore the average in Eq. (36) is equivalent to using a smeared source that has a smaller 
coupling to the opposite-parity excited states. This averaging procedure can be iterated in 
order to make the oscillating terms arbitrarily small. Empirically, we find that two iterations 
are sufficient for all of our numerical data: 
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(37) 



At our various light quark masses and lattice spacings the mass-splittings lie in the range 
0.1 ^ AmB ^ 0.3 in lattice units; thus use of the iterated average in Eq. (37) reduces 
the leading oscillating state amplitude by a factor of ~50-400 such that it can be safely 
neglected. 

In the case of the B ^ it 3-point correlation function, we wish to reduce both the 
oscillating contributions and the less visible non-oscillating contributions that arise from 
the cross-term between the lowest-lying pion and i?-meson opposite-parity states. If these 
contributions are reduced sufficiently, we can safely neglect all of them when extracting the 
ground-state amplitude A^^. We therefore construct a slightly more sophisticated average 
which combines the correlator both at consecutive time slices (t and t + 1) and at consecutive 
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source-sink separations (T and T + 1): 

-h^ t p-rrig (i -t) 

ct:\t,T) = — — 



X 



2C3^^7'^(t + l,T + l) C3^7'^(t + 2,r) C3^;^'^(t + 2,r+l 



g-E^O) g-m^,"' (T-t) g-^O) {t+2)^-mf {T-t-2) ^-E^°^ {t+2) ^-m%^ (T-t-l) 

+ 0{AEl AE^AniB, Ami). (38) 

This average reduces the unwanted parity states' contamination significantly. It ehminates 
both the leading 0{1) and sub leading (9(A£'^) contributions to the oscillating A^^ term, the 
two lowest-order (9(1, Am^) contributions to the oscillating A^^ term, and the 0{1, AEt^) 
contributions to the non-oscillating A^^ term. The size of the remaining A^^ term is a factor 
of ~7-20 times smaller than in the unsmeared 3-point correlator. 

We can now safely ignore contamination from opposite-parity states and determine the 
lattice form factors in a simple manner. We construct the following ratio of the smeared 
correlators: 



Cz7it,T) 2E, 



The lattice form factors are then: 



C;(t)C?(T-t)Ve-^^"^*e-™^^^-) 



= R;,o (t^T) (40) 



/r = ^RfrM. (41) 

We fit /||'^* and as defined in Eqs. (40)-(41) to a plateau in the region <C t ^ T where 
ordinary excited state contributions can be neglected. Figure 6 shows the determinations 
of /||'^* (left plot) and (right plot) for all of the momenta that we use in the chiral 
extrapolation on the coarse ensemble with ami /arris = 0.02/0.05. In practice, we fit a range 
of four time slices, choosing the interval that results in the best correlated confidence level. 
We have cross-checked the determination of the form factors via Eqs. (40)-(41) against 
determinations of the form factor that explicitly include excited state dependence in the fit 
ansatz and find that the results agree within errors. Our preferred method, however, yields 
the smaller statistical uncertainty in the form factors. 
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FIG. 6: Determination of the form factors /y (left plot) and f± (right plot) from plateau fits to 
the ratios defined in Eqs. (40) and (41). The statistical errors on the data points are from a single- 
elimination jackknife. The statistical errors in the plateau determination are from separate fits of 
500 bootstrap ensembles. 

C. Continuum and chiral extrapolation 

The quark masses in our numerical lattice simulations are heavier than the physical up 
and down quark masses. The effects of non-zero lattice spacings in Asqtad simulations are 
also too large to be neglected. In order to account for these facts, we calculate the desired 
hadronic matrix elements for multiple values of the light quark masses and lattice spacing, 
and then extrapolate to the physical quark masses and continuum using functional forms 
from heavy-light meson staggered chiral perturbation theory (HMSxPT) [42] . The HMSxPT 
expressions are derived using the symmetries of the staggered lattice theory, and therefore 
contain the correct dependence of the form factors on the quark mass and lattice spacing. 
In the case of the B niu form factors, the HMSxPT expressions are also functions of the 
pion energy (recall that we work in the frame where the S-meson is at rest). 

HMSxPT is a systematic expansion in inverse powers of the heavy quark mass. In the 
chiral and soft pion limits (m; — > and E,^ 0), the leading-order continuum HM^PT 
expressions for /y and f± take the following simple forms: 

hiE.) = ^ (42) 

J TV 
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where (j)B = fB\/'^B-, Jb is the S- meson decay constant, and fj^ is the pion decay constant. 
The coefficient Qb^Btt parameterizes the size of the B*-B-7i couphng. In the static heavy 
quark hmit, heavy quark spin symmetry does not distinguish between the pseudoscalar B- 
meson and the vector i?*-meson, which imphes that the decay constant (pB* = 4'b and the 
mass difference = tub* —ttlb 0. Inclusion of the parameter A^, however, ensures the 
proper location of the pole at m^, in the physical form factor At the next order in 

the heavy quark expansion, 0{l/mb) corrections split the degeneracy between the B- and 
5*-meson masses and decay constants. Furthermore, in the chiral and soft pion limits, all 
l/rrib corrections can be absorbed into the values of the parameters 0^, 0b* , Qb'B-k, and 
A^ [66]; thus /y and /j_ retain the functional forms in Eqs. (42) and (43) even at NLO in 
HMxPT. 

At lowest-order in S^PT, discretization effects split the degeneracies among the 16 tastes 
of pseudo-Goldstone mesons: 



"^L,s = i^x + ruy) + Aj 



(44) 



where x and y indicate the quark flavors, /i is a continuum low-energy constant, and A= is 
the mass-splitting of a meson with taste H. An exact U{1)a symmetry protects the taste 
pseudoscalar meson from receiving a mass-shift to all orders in S^PT, implying that Ap = 0. 
In addition, at O{o?), a residual 5*0(4) taste-symmetry preserves the degeneracies among 
mesons that are in the same irreducible representation: P,V, A,T, I [26]. Numerically, the 
size of the taste-splittings turn out to be comparable to those of the pion masses for the 
a = 0.09 fm and a = 0.12 fm Asqtad staggered lattices used in this work [16]. 

We extrapolate our numerical form factor data using HMSxPT expressions derived to 
zeroth order in l/m^,. The fit functions therefore depend upon the three remaining expansion 
parameters: mi, a, and E.,^. The HMSxPT expressions for the form factors to 0{mi, a^, E"^) 
are given explicitly in Eqs.(65)-(67) of Ref. [42]. Schematically, they read 
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1 

f±{mi,E^,a) = -j- 
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where "logs" indicate non-analytic functions of the pseudo-Goldstone meson masses, e.g., 
m^ln(m^/A^). The continuum low-energy constant Qb-'B-k enters these expressions in the 
coefficients of the chiral logarithms, which are completely fixed at this order. We use the 
phenomenological value of Qb'Btt = 0.51 [44] for the central value and vary Qb^Btt by a 
reasonable amount (see Sec. IV B) to estimate its contribution to the systematic uncertainty. 
Because the size of the mass-splitting A*^ is poorly determined from the lattice data and 
is consistent with the physical value within statistical errors, we fix A*^ to the PDG value, 
45.78 MeV [13], in our fits. The chiral logarithms also depend upon six extra constants that 
parameterize discretization effects due to the light staggered quarks: the four taste splittings 
a^Ay, a^A^, a'^A^, a?Ai and the two flavor-neutral "hairpin" coefficients a^S'y and o?5'j^ [27]. 
These parameters can be determined separately from flts to light pseudoscalar meson masses 
and decay constants; we therefore hold them flxed to the values determined in Ref. [67] while 
performing the continuum-chiral extrapolation. The variation of these parameters within 
their statistical errors results in a negligible change to the extrapolated form factors. The 
flve terms analytic in m/, a^, and E-,^ absorb the dependence upon the scale in the chiral 
logarithms, A^, such that the form factor is scale-independent. We leave the tree-level 
coefficients C||°]_ and the NLO analytic term coefficients C||^]_-C||^^ as free parameters to be 
determined via the fit to the lattice form factor data. In practice, we omit the analytic term 
proportional to {2mi + rris) from our fits because the strange sea quark mass is tuned to 
approximately the same value on each of our ensembles and we have simulated only full QCD 
points. This term is therefore largely indistinguishable from the analytic term proportional 
to mi. We have checked that omission of the sea quark mass analytic term has a negligible 
impact on the form factors in the chiral and continuum limits. 

In both earlier unquenched analyses of the B Triu semileptonic form factor [33, 34], 
the chiral extrapolation is performed as a two-step procedure: first interpolate the lattice 
data to fiducial values of E^^ and then extrapolate the results to the physical quark masses 
and continuum independently at each value of E.,^. The function used for the interpolation 
(which is different in the two analyses) introduces a systematic uncertainty that is difficult 
to estimate. In both cases, the chiral-continuum extrapolation makes use of the correct 
functional forms derived in HMSxPT, but, by extrapolating the results for each value of i?^ 
separately, the constraint that the low-energy constants of the chiral effective Lagrangian 
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are independent of the pion energy is lost. This omission of valuable information about 
the form factor shape introduces a further error that is unnecessary. The new analysis 
presented here instead employs a simultaneous fit using HMSxPT to our entire data set 
(all values of mi, a, and E^^) to extrapolate to physical quark masses and the continuum 
and interpolate in the pion energy [68]. This improved method eliminates the systematic 
uncertainties introduced in the two-step interpolate-then-extrapolate procedure, and exploits 
the available information in an optimal way. 

We perform our combined chiral and continuum extrapolation using the method of con- 
strained curve fitting [69]. Although we know that lattice data generated with sufficiently 
small quark masses and fine lattice spacings, and, in the case of the B Tiiu form factor, 
sufficiently low pion energies, must be described by lattice xPT, we do not know precisely 
the range of validity of the effective theory. Furthermore, the order in xPT to which we must 
work and the allowed parameter values depend upon both the quantity of interest and the 
size of the statistical errors. We therefore need a fitting procedure that both incorporates 
our general theoretical understanding of the suitable chiral effective theory and accounts for 
our limited knowledge of the values of the low-energy constants and sizes of the higher-order 
terms. Constrained curve fitting provides just such a method. 

Next-to-leading order xPT breaks down for pion energies around and above the kaon 
mass. Less than half of our numerical form factor data, however, is below this cutoff. 
Therefore, although we do not expect NLO HMSxPT to describe our data through momen- 
tum p = 27r(l, 1,0)/L, we cannot remove those points without losing the majority of our 
data. Nor can we abandon the NLO HMS^PT expressions for /y and f±, Eqs. (45) and (46), 
which are the only effective field theory guides that we have for extrapolating the numerical 
lattice form factor data to the continuum and physical quark masses. We therefore perform 
the continuum-chiral extrapolation using the full NLO HMSxPT expressions for /y and f±, 
including the 1-loop chiral logarithms, plus additional NNLO analytic terms to allow a good 
fit to the data through p = 2tt{1, 1,0)/ L. The NNLO terms smoothly interpolate between 
the region in which xPT is valid and the region in which the pion energies are too large and 
the higher-order chiral logarithms in E^^ can be approximated as polynomials. 

We express the analytic terms in the formulae for /y and f±, Eqs. (45) and (46), as 
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products of dimensionless expansion parameters: 

0.05-0.19 (47) 
0.03-0.09 (48) 

0.22-0.78, (49) 

where A is the average staggered taste-sphtting and we show the range of values for each 
of these parameters corresponding to our numerical lattice data. (Note that we omit the 
p = 27r(l, 1, 1)/L data points from our chiral fits because these would lead to ~ 1-) 
Because each of the above expressions is normalized by the chiral scale, ^ 4:'^fn, the 
undetermined coefficients C||^]_-C||^]_ should be of 0{1) in these units. We therefore constrain 
the values of the low-energy constants C||''^-C||^^ in our fits with Gaussian priors of width 2 
centered about 0. 

The statistical errors in the numerical lattice data come from the 3-point fits described 
in the previous subsection. In order to account for the correlations among the various pion 
energies on the same sea quark ensemble in the chiral-continuum extrapolation, we preserve 
the bootstrap distributions. We perform a separate correlated fit to each of the 500 bootstrap 
ensembles in which we remake the full bootstrap covariance matrix for each fit. We average 
the 68% upper and lower bounds on the form factor distributions to determine the statistical 
and systematic errors in /y and f± that are plotted in Fig. 7 and presented in Table II below. 

Because we do not know a priori how many terms are necessary to describe the available 
lattice data, we begin with strictly NLO fits using the formulae for /y and f± in Eqs. (45) 
and (46). We fit the lattice data for /y and f± separately even though the ratio of leading- 
order coefficients, c^±^ /c^^\ is predicted to equal gB*BTT to NLO in xPT; this is because the 
value of Qb* b-k is known to only ~ 50% from phenomenology. We obtain a good fit of the /j_ 
lattice data to the NLO expression without the inclusion of higher-order NNLO terms. This 
is probably because the shape of is dominated by the 1 / [E^^ + /\*) behavior and therefore 
largely insensitive to the other terms. We cannot, however, obtain a good fit of f\\ to the 
strictly NLO expression, and must add higher-order terms in order to obtain a successful 
fit. Specifically, NNLO analytic terms proportional to miEj^ and are both necessary to 
achieve a confidence level better than 10%. 

Although we could, at this point, choose to truncate the HMS^PT extrapolation formulae 
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to include only those terms necessary for a good confidence level, we instead include "extra" 
NNLO analytic terms to both the /y and f± fits, constraining the values of their coefficients 
with Gaussian priors of ± 2. The introduction of more free parameters increases the 
statistical errors in the extrapolated values of the form factors; these larger errors reflect the 
uncertainty in the size of the newly-included higher-order contributions. We continue to add 
higher-order analytic terms until the central values of the extrapolated form factors stabilize 
and the statistical errors in the form factors reach a maximum. This indicates that any 
further terms are of sufficiently high order that they do not affect the fit and can safely be 
neglected. We find that this occurs once the extrapolation formulae for /y and f± contain all 
eight sea-quark mass-independent NNLO analytic terms. The inclusion of NNNLO analytic 
terms does not further increase the size of the error bars. 

Figure 7 shows the preferred constrained fits of f\\ (upper plot) and f± (lower plot) 
versus E^, where both the x- and y-axes are in ri units. ^ Each fit is to the NLO HMSxPT 
expression, Eqs. (45) and (46), plus all sea-quark mass-independent NNLO analytic terms. 
The square symbols indicate fine lattice data, while the circles denote coarse data. The 
six colored curves show the fit result projected onto the masses and lattice spacings of the 
six sea quark ensembles; the red line should go through the red circles, and so forth. The 
thick black curve shows the form factor in the continuum at physical quark masses with 
symmetrized bootstrap statistical errors. 

We use functions and constraints based on HMSxPT to perform the chiral-continuum 
extrapolation because we know that HMSxPT is the correct low-energy effective description 
of the lattice theory. Nevertheless, we must compare various properties with theoretical 
expectations in order to check for overall consistency. An essential first test is that we 
can successfully fit the data with good confidence levels and obtain low-energy coefficients 
that are of the predicted size. We can also verify the convergence of the series expansion 
by calculating the ratios of the higher-order contributions to the leading-order form factor 

^ As a cross-check of the constrained fits, we also perform unconstrained fits of /|| and f± with only the 
minimal number of analytic terms needed for a good fit. The results are consistent, but the unconstrained 
fit results have smaller statistical errors because they include 6-8 fewer fit parameters. 
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xVd.o.f. = 1.1 
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ami /anis 
ami /ams 
ami jams 
ami lams 
ami jams 
ami jams 



— continuum QCD 



0.0062/0.031 fine 
0.0124/0.031 fine 
0.005/0.05 coarse 
0.007/0.05 coarse 
0.01/0.05 coarse 
0.02/0.05 coarse 



1.5 




xVd.o.f. = 0.32 
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FIG. 7: Chiral-continuum extrapolation of /y (upper) and /_l (lower) using constrained NLO 
HMSxPT plus all NNLO analytic terms with gB'B-w = 0.51 and ri = 0.311 fm. The square 
symbols indicate a ~ 0.09 fm lattice data points while the circular symbols indicate a ~ 0.12 fm 
coarse data points. The black curve is the chiral-continuum extrapolated form factor symmetrized 
bootstrap statistical errors only. 
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contributions: 



/l 
/i 



NLO 



NNLO 



fLO 

=500 MoV J ± 



/l 



LO 



47%, 
3%, 



NLO 



/: 

jNNLO 



=500 MeV 



LO 



=500 MeV 



=500 McV 



48%, (50) 
4%, (51) 



where we choose a nominal value of E^^ = 500 MeV for illustration because it is on the high 
end of the expected range of validity of xPT. Finally, because the leading-order coefficient, 
C||°'', is expected to be equal to = fs^Jf^B in HM^PT, we can compare its value with 
that oi ipB determined from our preliminary decay constant analysis. Although the B- 
meson decay constant calculation uses the same staggered gauge configurations, it employs 
different heavy-light meson 2-point correlation functions with the axial current, a different 
HMSxPT fit function, and different perturbative renormalization factors, and is therefore 
largely independent of the B niu semileptonic form factor calculation. For the preferred 
/|l fit shown in Fig. 7, we find c^^^ = 0.81 ± 0.07, where the errors are statistical only. 
This is quite close to our current preliminary result, r'^^'^cpB = 0.92 ± 0.03 (statistical error 
only) [70, 71], especially considering that the HMSxPT extrapolation formula for f\\ neglects 
some of the 0{l/mij) contributions. 

An interesting use of our numerical B — > ttIv form factor data is to determine the 
approximate value of the B*-B-'k coupling, Qe'Ett, from lattice QCD. For the preferred fits 
shown in Fig. 7, we find that the ratio of leading-order coefficients is 

gB'B. ~ = 0.22 ± 0.07, (52) 

C|| 

and is independent of the choice for Qb'Btt in the chiral logarithms within statistical errors. 
This determination omits the 0{l/mb) corrections to the chiral logarithms in the HMSxPT 
extrapolation formulae for /y and f±, Eqs. (45) and (46), and neglects the difference between 
(pB and (pB*] we do not attempt to estimate the systematic uncertainty introduced by these 
or other effects. The value is lower than the determination of Stewart, gB*B-K = 0.51, which 
comes from a combined analysis of several experimental quantities, including the D*-meson 
decay width, through (9(l/mc) in HM^PT [44]. It is consistent, however, with the range of 
values determined by the HPQCD Collaboration, who allowed Qb'Btt to be a free parameter 
in their chiral-continuum extrapolation and found < gB*Bn ^ 0.45 [33]. 
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IV. ESTIMATION OF SYSTEMATIC ERRORS 



In this section, we discuss all of the sources of systematic uncertainty in our calculation of 
the B -H> 7t£u form factor fj^{q'^). We present each error in a separate subsection for clarity. 
The value of the form factor along with the total error budget, is given in Table II. 

A. Chiral-continuum extrapolation fit ansatz 

We use the method of constrained curve fitting to estimate the effect of neglected higher- 
order terms in the HMSxPT chiral-continuum extrapolation formulae. Our fit procedure is 
described in detail in Sec. Ill C. Therefore, the errors in f\\ and extrapolated to physical 
quark masses and zero lattice spacing shown in Fig. 7 reflect both the statistical errors 
in the Monte Carlo data and the systematic errors due to our limited knowledge of the 
higher-order terms, which we specified with priors. We do not need to include a separate 
systematic uncertainty due to the choice of fit function, as would be the case had we used 
an unconstrained fit with fewer terms. 

B. Qb'Btt uncertainty 

We fix the size of the B*-B-tx coupling to Qb'Btt = 0.51 in the coefficients of the chi- 
ral logarithms while extrapolating to the physical light quark masses and continuum. Our 
choice is based upon the following considerations. Because the coupling gB'Bn is expected 
to be approximately equal to gD*Dn due to heavy-quark symmetry, we use the phenomeno- 
logical value of the D*-D-ti coupling determined by Stewart in Ref. [44], which comes from 
a combined analysis of several experimental quantities that includes the D*-meson decay 
width [72]. This value is presented without errors, and is an update of Stewart's earlier 
analysis in Ref. [73] which incorporates additional experimental results. His earlier calcu- 
lation finds a significantly lower value of go'Dw = 0-27^qq2-o where the first errors are 
experimental and the second errors come from an estimate of the sizes of the 1-loop coun- 
terterms [73]. A more recent phenomenological determination of the D*-D-tt coupling by 
Kamenik and Fajfer, which also includes up-to-date experimental data, improves upon the 
analysis method of Stewart by including contributions from both positive and negative parity 
heavy mesons in the loops [74]. They find an even higher value of Qd'O-k = CGG^^^ggg, where 
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the uncertainty only reflects the error due to counterterms. We therefore conclude that, al- 
though recent experimental measurements of the D* width may constrain the coupling [72] 
at tree-level, the size of Qb'B-k is not well-determined in the literature. 

In order to determine the error in the form factor from the uncertainty gB*B-K we vary 
the parameter over a generous range. The smallest value of qb* b-k that we have seen in the 
literature is go'Dw = 0.27 [73]. The largest is go'Dn = 0.67, which comes from a quenched 
lattice calculation [75]. (There has not yet been an unquenched "2+1" flavor determination 
of Ob* Btt-) We therefore vary qb* b-k over the entire range from 0.27-0.67 and take the largest 
difference from the preferred determination of f+{q^) using gB*Bn = 0.51 as the systematic 
error due to the uncertainty in the B*-B-it coupling. The lattice data is largely insensitive 
to the value of gB*Bn in the coefficient of the chiral logarithms; all values of the parameter 
yield similar fit confldence levels. The resulting systematic uncertainty in f+{q'^) is less than 
3% for all bins despite varying gB*Bn by almost 50%. 

C. Scale (ri) uncertainty 

We use the MILC Collaboration's determination of the scale from their calculation of /tt, 
Ti = 0.311 fm, to convert between lattice and physical units [52]. The parameter ri enters 
the form factor calculation in a number of places: we use the PDG values of and in the 
chiral- continuum extrapolation formulae [13], we set m^r to the PDG value in the resulting 
flt functions to determine the form factors at the physical point, and we convert /y and fj_ to 
physical units via ri before combining them to extract An alternative determination 

using the HPQCD Collaboration's lattice data for the T 2S'-1S' [76] splitting yields a result 
that is ~ 2% larger, ri = 0.317 fm. We therefore repeat the chiral-continuum extrapolation 
of /|| and f± using this higher value of ri, combine them into the dimensionless form factor 
using this higher value of ri, and take the difference from the preferred form factor 
result as the systematic error due to uncertainty in the overall lattice scale. The difference 
ranges from 1-1.5% for most values. This is consistent with our naive expectation that 
a ~ 2% difference in ri will result in a ~ 1% difference in f+{q^) because f\\ has dimensions 
of GeV^^^ and f±_ has dimensions of GeV"^''^. 
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D. Light quark mass (m, rUg) determinations 

We obtain the form factors f\\ and f± in continuum QCD by setting the lattice spacing to 
zero and the light quark masses to their physical values in the HMSxPT expressions, once 
the coefficients have been determined from fits to the numerical lattice data. We use the 
most recent calculations of the bare quark masses by the MILC Collaboration from fits to 
light pseudoscalar meson masses: 



where m is the average of the up and down quark masses and the quoted errors include both 
statistics and systematics [52]. We vary the bare light quark mass, rim, within its stated 
uncertainty and take the maximal difference from the preferred form factor result to be the 
systematic error; we find that the error is 0.3% or less for all values of q^. We perform the 
same procedure for the bare strange quark mass, and find that the resulting error ranges 
from ~ 0.5-1.5% over the various bins. 

E. Bottom quark mass (mt,) determination 

The value of the form factor f+{q^) depends upon the 6-quark mass, which we fix to its 
physical value throughout the calculation. Specifically, we first determine the value of the 
hopping parameter, k, in the SW action for which the lattice kinetic mass agrees with the 
experimentally-measured i^^-meson mass. We then use this tuned ki, when calculating all 
of the 2- and 3-point heavy-light correlators needed for the B nlv form factor. With our 
current tuning procedure we are able to determine to ~ 6% accuracy. This uncertainty 
in Kb is conservative; it is primarily due to poor statistics, and will decrease considerably 
after the analysis of the larger data set that is currently being generated. 

The uncertainty in ki, produces an uncertainty in the form factor. We estimate this by 
calculating the form factor /+ at two additional values of n (one above and one below the 
tuned value) on the ami/anis = 0.02/0.05 coarse ensemble. This is sufficient because the 
heavy-quark mass-dependence of the form factor is largely independent of the sea quark 
masses and lattice spacing. We find the largest dependence upon Hf, at momentum p = 



nm X 10^ = 3.78(16) 
rim, X 10^ = 102(4), 



(53) 
(54) 
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FIG. 8: Normalized form factor at momentum p = 27r(l,l,0)/L as a function of k on the 
anil/ arris = 0.02/0.05 coarse ensemble. The central data point corresponds to the tuned Hb, and 
the thick red line shows a linear fit to the three data points. The two dashed vertical lines indicate 
the upper and lower bounds on ki,. 

27r(l, 1, 0)/L, shown in Fig. 8, for which a 6% uncertainty in Kh produces a 1.2% uncertainty 
in the form factor. We therefore take 1.2% to be the systematic error in f+{q^) due to 
uncertainty in the determination of the 6-quark mass. 



F. Gluon and light-quark discretization errors 

We estimate the size of discretization errors in the form factor f+{(f) with power-counting. 
We choose conservative values for the parameters that enter the estimates: A = 700 MeV 
and av{q*) = 1/3, which is a typical value on the fine lattice spacing [58, 63]. 

We calculate the B txIv semileptonic form factor using a one-loop Symanzik-improved 
gauge action for the gluons [77-80] and the Asqtad-improved staggered action for the light 
up, down, and strange quarks [81, 82]. Because both the gluon and light quark actions are 
(9(a^)-improved, the leading discretization effects are of (9(as(aA)^). We parameterize these 
errors in the fit to numerical lattice form factor data by including analytic terms proportional 
to a? in the HMSxPT extrapolation formulae for /y and f±. Because we only have data at 
two lattice spacings, however, we do not include a separate term proportional to as{ah)'^ to 
account for the fact that differs by a few percent between the lattice spacings. We then 
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remove the majority of light quark and gluon discretization effects from the final result by 
taking a — > 0. Similarly, we identify and remove higher-order discretization effects in the 
chiral-continuum extrapolation through the NNLO analytic terms in the fit functions. The 
remaining gluon and light quark discretization errors are negligible. 

The calculation of the B iriu form factor requires 2-point and 3-point functions with 
nonzero momenta, which introduces momentum-dependent discretization errors. The lead- 
ing p-dependent discretization error is of Olasiap)"^)- We parameterize these errors, up to 
variations in at the two lattice spacings, with the two NNLO analytic terms proportional 
to a^E'^ and a'^rrii in the extrapolation formulae for /y and f± and remove them by taking the 
continuum limit of the resulting fit functions. This also largely removes errors of (9(a^(ap)^). 
The remaining momentum-dependent discretization effects are of 0{as{ap)^). On the 28^ 
fine lattices, as{ap)'^ = 0.003 for our highest-momentum data points with ap = 27r(l, 1, 0)/28. 
Therefore the uncertainty in f+{q^) due to momentum-dependent discretization effects is 
negligible compared with our other systematic errors. 



G. Heavy-quark discretization errors 

We use HQET as a theory of cutoff effects to estimate the size of discretization errors 
due to use of the Fermilab action for the heavy bottom quark. Because both the lattice 
and continuum theories can be described by HQET, heavy-quark discretization effects can 
be classified as a short-distance mismatch of higher- dimension operators [37-39]. Each 
contribution to the error is given by [83] 

error, = | [C^^^Q, rn,a) - Cr^niQ)] {O,) \ , (55) 

where Oi is an effective operator, and Cj'^^ {mq , moa) and Cf°'^*(mQ) are the corresponding 
short-distance coefficients when HQET is used to describe lattice gauge theory or continuum 
QCD, respectively. The coefficient mismatch can be written as 

C]^\mQ, moa) - Cr'img) = a'^'-^-V.(rnoa), (56) 

and the relative error in our matrix elements can be estimated by setting (Oi) ~ A^^^''^. 
Then each contribution to the error is 

error, = /,(moa)(aAQCD)"''" (57) 
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recovering the counting in powers of a familiar from Symanzik, while maintaining the full 
moa dependence. The functions fi can be deduced from Refs. [36, 84] and are compiled 
in Appendix A. Adding all contributions of 0{asa) and O(a^) from the action and the 
current, we obtain a relative error of 2.84% (4.16%) for /y and 3.40% (4.98%) for fj_ on the 
fine (coarse) lattices. We therefore take 3.4% to be the error in f+{q^) due to heavy-quark 
discretization effects. 

H. Heavy-light current renormalization 

We determine the majority of the heavy-light current renormalization nonperturbatively. 
The dependence of Zy on the sea quark masses and on the mass of the light spectator 
quark in the 3-point correlator are both negligible; the statistical error in Zy is ~ 1%. The 
dependence of Zy on the sea quark masses is also negligible, and the statistical error in Zy 
is ~ 1%. We therefore include a/ (1%^ + 1%^) ~ 1.4% as the systematic uncertainty in the 
form factor f+{q^) due to the uncertainty in the nonperturbative renormalization factors Zy 
and Zy for all values of g^. 

We determine the remaining renormalization of the heavy-light current using lattice per- 
turbation theory. The 1-loop correction to fj_ is ~ 3% on the fine ensembles and ~ 4% 
on the coarse ensembles. Because we calculate py^ to 0{as), the leading corrections are of 
0{a1). We might therefore expect the 2-loop corrections to py^ to be a factor of a.^ smaller, 
or ~ 1%. In order to be conservative, however, we take the entire size of the 1-loop correc- 
tion on the fine lattices, or 3%, as the systematic uncertainty in f+{q^) due to higher-order 
perturbative contributions for all bins. 

I. Tadpole parameter (mq) tuning 

In order to improve the convergence of lattice perturbation theory, we use tadpole- 
improved actions for the gluons, light quarks, and heavy quarks [48]. We take Uq from 
the average plaquette for the gluon and sea quark action [18]. On the fine lattice, we make 
the same choice for the valence quarks. For historical reasons, however, we use Uq determined 
from the average link in Landau gauge for the valence quarks on the coarse ensembles. The 
difference between uq from the two methods is 3-4% on the coarse ensembles. We must. 
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therefore, estimate the error in the form factor due to this poor choice of tuning. 

The tadpole-improvement factor enters the calculation of f+{q^) in several ways. The 
factor of Mo that enters the normalization of the heavy Wilson and light staggered quark 
fields cancels exactly between the {ti\V^\B) lattice matrix element and the nonperturbative 
renormalization factor a/ ZyZy. The most significant effect of the mixed uq values is in the 
chiral-continuum extrapolation of /y and f±. The different choices for valence and sea quark 
actions imply that the coarse lattice data is partially quenched. We study this effect by 
performing the chiral extrapolation in two ways: one assuming that both valence and sea 
quarks have the mass of the sea quark and the other assuming that both have the mass of 
the valence quark. This leads to a 3% error in the highest bin, and a ~ 1 — 1.5% error in 
the bins that affect the determination of \Vub\- Most of the other effects of changing uq in the 
lattice action and current can be absorbed into our estimate of the uncertainty from higher- 
order perturbative corrections to py, to discretization errors, and to the normalization of 
the Naik term. All but the last are already budgeted in Table II. The Naik term in the 
Asqtad action ensures that the leading discretization errors in the pion dispersion relation 
are 0{asa'^p'^)- We therefore estimate the error in f+{q'^) due to different Naik terms in the 
valence and sea sectors to be equal to the largest value of agO^p^ on the coarse lattice times 
the ratio of the Landau link over plaquette Uq cubed, or ~ 0.2%. 

We add the flat error from the Naik term to the bin-by-bin error due to the light quark 
mass used in the chiral extrapolation in quadrature to obtain the total uncertainty. Although 
our estimate is of necessity rather rough, we find that the errors due to uq tuning are much 
smaller than the dominant errors in Our error estimate is therefore adequate for the 

determinations of the B iriu form factor and \Vub\ presented in this work. 

J. Finite volume effects 

We estimate the uncertainty in the form factor f+{q'^) due to finite- volume effects using 
1-loop finite volume HMSxPT. The finite volume corrections to the HMSxPT expressions 
for /|| and are given in Ref . [42] in terms of integrals calculated in Ref . [85] . It is therefore 
straightforward to find the relevant corrections for our simulation parameters. We find that 
the 1-loop finite volume corrections are well below a percent for all of our lattice data 
points. Because finite volume errors increase as the light quark mass decreases, they are 
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largest on the ami jams = 0.007/0.05 coarse ensemble. The biggest correction is to f± at 
p = 27r(l, 1,0)/L, and is 0.5%. We therefore take this to be the uncertainty in f+{q^) due 
to finite volume errors for all bins. 



V. MODEL-INDEPENDENT DETERMINATION OF iKbl 

It is well-established that analyticity, crossing symmetry, and unitarity largely constrain 
the possible shapes of semileptonic form factors [86-89] . In this section we apply constraints 
based on these general properties to our lattice result for the form factor f+iq"^) and thereby 
extract a model-independent value for the CKM matrix element \Vub\- 

Until now the standard procedure used to extract |\4fe| from B niu semileptonic decays 
has been to integrate the form factor over a region of g^, and then combine the 

result with the experimentally measured decay rate in this region: 

w ^ C "'^ K'"^ ^ - '^'^ - ""^^--^''^ '''' 

^ ^min 

The integration, however, necessitates a continuous parameterization of the form factor over 
the full range from g^j^ to g^^^. 

In our earlier, preliminary unquenched analysis, we determine /+(g^) by fitting the lattice 
data points to the Becirevic-Kaidalov (BK) parameterization [90], 



Uil') = n (59) 



_ /+(0) 

[1 — (p) {1 — a(p 

2> _ /+(0) 



where q^ = g^/m^,. The BK ansatz contains three free parameters and incorporates many 
of the known properties of the form factor such as the kinematic constraint at g^ = 0, 
heavy-quark scaling, and the location of the B* pole. The HPQCD Collaboration instead 
uses the four-parameter Ball-Zwicky (BZ) parameterization [91], which is the same as the 
BK function in Eq. (59) plus an additional pole to capture the effects of multiparticle states. 
In both cases, however, the choice of fit function introduces a systematic uncertainty that 
is difficult to quantify. 

It is likely the BK and BZ parameterizations can be safely used to interpolate between 
data points, whether they be at high g^ from lattice QCD or at low g^ from experiment. 
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TABLE II: Statistical and systematic error contributions to the B ttIu form factor. Each source of uncertainty is discussed in Sec. IV. 
For each of the 12 bins, the error is shown as a percentage of the total form factor, which is given in the second row from the 

top. Because the bootstrap errors in the form factor are asymmetric, the errors shown are the average of the upper and lower bootstrap 
errors. In order to facilitate the use of our result, we also present the normalized statistical and systematic bootstrap correlation matrices 
in Table IV and the total bootstrap covariance matrix in Table V. 
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It is less clear, however, how well these ansatze can be trusted to extrapolate the form 
factor shape beyond the reach of the data points. Furthermore, comparisons of lattice and 
experimental determinations of BK or BZ fit parameters are not necessarily meaningful. 
For example, if the slope parameters a from experiment and lattice QCD were found to 
be inconsistent, we would not know whether theory and experiment disagree, or whether 
the parameterization is simply inadequate. A parameterization that circumvents this issue 
is therefore desirable. In this work we pursue an analysis based on the model-independent 
2;-parameterization, which is pedagogically reviewed in Ref. [87]. 



A. Analyticity, unitarity, and heavy-quark constraints on heavy-light form factors 

All form factors are analytic functions of (f except at physical poles and threshold branch 
points. In the case of the B irlu form factors, fiq"^) is analytic below the Btt production 
region except at the location of the B* pole. The fact that analytic functions can always be 
expressed as convergent power series allows the form factors to be written in a particularly 
useful manner. 

Consider mapping the variable onto a new variable, z, in the following way: 

z{,\ to) = V^^Z.V^^ (61) 

v^-<iyu + v^-to/u 

where = (m^ + m^)^, t_ = [itlb — rrij^y, and to is a free parameter. Although this 
mapping appears complicated, it actually has a simple interpretation in terms of g^; this 
transformation maps > t+ (the production region) onto |2| = 1 and maps < t+ (which 
includes the semileptonic region) onto real z G [— 1, 1]. In terms of z, the form factors have 
a simple form: 

where the Blaschke factor P{q^) is a function that contains subthreshold poles and the outer 
function 0(5^, to) is an arbitrary analytic function (outside the cut from < < oo) whose 
choice only affects the particular values of the series coefficients a^. 

For the case of the B iriu form factor f+^q"^), the Blaschke factor P+(g^) = z{q^^m\,) 
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accounts for the B* pole. In this work we use the same outer function as in Ref. [44]: 

3/2 

Mi'^to) = ^|—^(^/t+-q' + Vt+-to]{^/U-q' + Vt+-t ' 




(63) 



where x^j^ is a numerical factor that can be calculated via the operator product expansion 
(OPE) [88, 92]. This choice of 0+(g^,to); when combined with unitarity and crossing- 
symmetry, leads to a particularly simple constraint on the series coefficients in Eq. (62). 
Although the t-dependence of Eq. (63) appears complicated, it is designed so that the sum 
over the squares of the series coefficients is t-independent: 

Y.^l = l^A i\nz)<P{z)f\z)\' ^ A, (64) 

fc=0 

where the value of the constant A depends upon the choice of Xj''' in Eq. (63). Because the 
decay process B — » ttIv is related to the scattering process Btt by crossing symmetry, 
the sum of the series coefficients is bounded by unitarity, i.e., the fact that the production 
rate of Bti states is less than or equal to the production of all final states that couple to the 
h ^ u vector current. In particular, if one chooses the numerical factor Xj"* to be equal to 
the appropriate integral of the inclusive rate — > X;,, the sum of the coefficients is bounded 
by unity: 

TV 

Y.al<l, (65) 

fc=0 

where this constraint holds for any value of N and the symbol indicates higher-order 
corrections to x^j^ iii ots and the OPE. 

Such higher-order corrections turn out to be negligible for the B niu form factor 
because the bound in Eq. (65) is far from saturated, i.e., the sizes of the coefficients turn 
out to be much less than one. Becher and Hill [93] have pointed out that this is due to the 
fact that the 6-quark mass is so large. In the heavy-quark limit, the leading contributions 
to the integral in Eq. (64) are of (9(A^/m^), where A is a typical hadronic scale. Assuming 
that the ratio A/rrib ~ 0.1, the heavy-quark bound on the a^s is approximately thirty times 
more constraining than the bound from unitarity alone: 



-O-OOl- (66) 



fc=0 
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We point out that the authors of Ref. [46] have recently proposed a shghtly different 
parameterization of the B niu form factor with a simpler choice of outer function, 0=1: 

^ oo 

ui^') = r— v;;^E^^(^o)^(^''^o)^ (67) 

This choice enforces the correct scaling behavior, f+iq"^) ~ as ^ oo. It leads, 

however, to a more complicated constraint on the series coefficients: 

N 

J2 Bjkhjhk < 1, (68) 

i,fc=0 

where the elements of the symmetric matrix Bj^ are calculable functions of t^. Because 
B iriu semileptonic decay is far from — > oo, and because the unitary bound is so far 
from being saturated, the choice of outer function should make a negligible impact on the 
resulting determination of \Vub\- We therefore use the more standard outer function given 
in Eq. (63) because the constraint in Eq. (65) is independent of the number of terms in the 
power series, and is therefore simpler to implement. 

The free parameter to can be chosen to make the maximum value of as small as possible 
in the semileptonic region; we choose to = 0.65t_ as in Ref. [44]. For B —* ttIu semileptonic 
decays this maps the physical region onto: 

< t < t_ ^ -0.34 < z < 0.22. (69) 

The bound on the coefficients in the z-expansion combined with the small numerical values 
of \z\ in the physical region ensures that one needs only the first few terms in the ^-expansion 
to accurately describe the form factor shape. Moreover, as the precision of both the lattice 
calculations and experimental measurements improve, one may easily include higher-order 
terms as needed. 

B. Determination of \Vub\ using z-parameterization 

In 2007 the BABAR Collaboration published a measurement of the shape of the B — > iriu 
semileptonic form factor with results for 12 separate bins between g^j^ ^ 1 GeV^ and 
?max ~ 24 GeV^ [43]. This suggests that lattice QCD calculations are now needed primarily 
to provide a precise form factor normalization at one value of g^ in order to determine \Vub\- 
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The minimal error in \Vub\ can, of course, still be attained by using all of the available 
information on the form factor shape and normalization, provided that one analyzes the 
data in a model-independent way. 

Because as many terms can be added as are needed to describe the B niu form factor 
to the desired accuracy, use of the convergent series expansion allows for a systematically 
improvable determination of \Vub\- We fit our lattice numerical Monte Carlo data and the 
12-bin BABAR experimental data together to the 2;-expansion, leaving the relative nor- 
malization factor, \Vub\, as a free parameter to be determined by the fit. In this way we 
determine \Vub\ in an optimal, model-independent way. 

We first fit the lattice numerical Monte Carlo data and the 12-bin BABAR experimental 
data separately to the z-expansion in order to check for consistency. We use Gaussian 
priors with central value and width 1 on each coefficient in the ^-expansion to impose the 
unitarity constraint. Although this manner of constraining the coefficients is less stringent 
than the strict bound given in Eq. (65), the choice does not matter because the unitary 
bound is far from saturated and the individual coefficients all turn out to be much less than 
1. We obtain identical fit results even when the coefficients are completely unconstrained. 

The left-hand plot in Fig. 9 shows the BABAR measurement of the B niu semileptonic 
form factor, f+{(f) [43]. The right-hand plot shows the same data multiplied by the functions 
P+(g^) and (j)^{q'^,tQ) and plotted versus the variable z. After remapping from to z 
there is almost no curvature in the experimental data. This indicates that most of the 
curvature in the data is due to well-understood QCD effects that are parameterized by 
the functions P+(g^) and 0+(g^,to)- Consequently the experimental data is well- described 
by a normalization (oq) and slope (ai/ag), as shown in Fig. 9. The slope of the BABAR 
experimental B — > iriu form factor data is 

— = -1.60 ±0.26. (70) 
ao 

If one includes a curvature term in the z-fit, the coefficient a2 is poorly determined, but is 
found to be negative at ~ 1.5a. The value of Oi is consistent with the result of the linear 
fit. 

Figure 10 shows the lattice determination of the B iriu semileptonic form factor, /_|_ 
vs. (left plot) and the remapped form factor, P^(f)^f\ vs. z (right plot). As is the case 
for the experimental data, the shape of the lattice form factor is less striking after taking 
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FIG. 9: Experimental data for the B ttIv form factor times the CKM element \Vub\ from the 
BABAR collaboration [43]. The left plot shows \Vub\ x /+ versus while the right plot shows 
\Vub\ X f+ multiplied by the functions P+(j)+ and plotted against the new variable z. Both the 
2-parameter fit (dashed blue line) and 3-parameter fit (solid red curve) have good x^/d.o.f.'s. 

out the B* pole and other known QCD effects. When the lattice calculation of the form 
factor is fit to the ^-parameterization, however, it determines both a slope and a curvature. 
One cannot, in fact, successfully fit the lattice data without including a curvature term. The 
slope and curvature of the lattice determination of the B iriu form factor are 



-1.75 ±0.91, 



ai 
ao 

— = -5.22 ±1.39. 

ao 



(71) 
(72) 



The above uncertainties are the standard errors computed from the inverse of the parameter 
Hessian matrix that result from a fit using the full covariance matrix determined from 
the bootstrap distributions of chiral-continuum extrapolated values of f\\ and f±, including 
systematics. 

Because the shapes of the lattice calculation and experimental measurement of the form 
factor are consistent, we now proceed to fit them simultaneously to the ^-expansion and 
determine \ Vub\- The numerical lattice and measured experimental data are independent, so 
we construct a block-diagonal covariance matrix where one block is the total lattice error 
matrix and the other is the total experimental error matrix. The combined fit function 
includes the series coefficients (a^'s) plus an additional parameter for the relative normal- 
ization between the lattice and experimental results (|V„fe|)- In order to account for the 
systematic uncertainty in \Vub\ due to poorly-constrained higher-order terms in z, we con- 
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FIG. 10: Lattice calculation of the B — > ttIu form factor. The left plot shows /+ vs. q'^ while the 
right plot shows P-^(j)^f^ vs. z. The inner error bars indicate the statistical error, while the outer 
error bars indicate the sum of the statistical and systematic added in quadrature. A 3-parameter 
z-fit is needed to describe the lattice data with a good x^/d.o.f. 

tinue to add terms in the series until the error in \ Vub\ reaches a maximum. This occurs once 
we include the term proportional to z^. The resulting combined 2- fit is shown in Fig. 11, 
and the corresponding fit parameters are 



V^b\ X lO'^ = 3.38 ±0.36, 

ao = 0.0218 ±0.0021, 

ai = -0.0301 ±0.0063^ 

02 = -0.059 ±0.032, 

as = 0.079 ±0.068. 



(73) 
(74) 
(75) 
(76) 
(77) 



The values of the coefficients are all much smaller than one, as expected from heavy-quark 
power-counting. The sum of the squares of the coefficients is Y2(^k ~ 0.011 ± 0.012, and is 
consistent with the prediction of Becher and Hill within uncertainties in the series coefficients 
and in the choice of the hadronic scale in Eq. (66) [93]. 

By combining all of the available numerical lattice Monte Carlo data and 12-bin BABAR 
experimental data for the B — ^ Triu form factor in a simultaneous fit we are able to determine 
\Vub\ to ~ 11% accuracy. This error is independent (within < 0.5%) of the choice of the 
parameter to used in the change of variables from to 2(g^,to) and in the outer function 
0+(g^,to)- In order to demonstrate the advantage of the combined fit method, we compare 
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xVd-oi- = 0.59 xVd-oi. = 0.59 




FIG. 11: Model-independent determination of \Vub\ from a simultaneous fit of lattice and exper- 
imental B nil' semileptonic form factor data to the z-parameterizaton. The left plot shows 
P^4>jf.f\ vs. z while the right plot shows /+ vs. q^. Inclusion of terms in the power-series through 
yields the maximum uncertainty in \ Vub\] the corresponding 4-parameter z-fit is given by the red 
curve in both plots. The circles denote the Fermilab-MILC lattice data, while the stars indicate the 
12-bin BABAR experimental data, rescaled by the value of \Vub\ determined in the simultaneous 
z-fit. 

the error in \Vub\ given in Eq. (73) with that obtained from separate z-fits of the lattice 
and experimental data. A 2;-fit to the 12-bin BABAR experimental data alone determines 
the normalization ao°^^ to ~ 8%, while a z-fA to our numerical lattice data determines 
a^^^ to ~ 14%. Thus separate fits lead to a determination of \Vub\ = cto'^^^/'^o^^* with an 
approximately 16% total uncertainty.^ The combined fit yields a significantly smaller error 
and is thus preferred. 

When the numerical lattice data and experimental data are fit simultaneously, utilizing 
all of the available data points is of secondary importance for reducing the total uncertainty 

^ Because the values of the coefficients of the power-series in z depend upon the choice of the parameter to in 
Eqs. (61)-(63), we could, in principle, choose a different value of to in order to minimize the error in either 
00*"^^ or ao^'^*. For example, use of to = 22.8 GeV^ reduces the uncertainty in the lattice normalization 
because the error in the lattice form factor is smallest at this g^-value. Use of to ~ 22.8 GeV^ greatly 
increases the uncertainty in the experimental normalization, however, because the experimental data is 
poorly-determined at large values of . Ultimately, this choice of leads to an even worse determination 
of iKibl than from our standard choice of — 0.65t_. Although we did not attempt to determine the 
value of <o that minimizes the total error in | V^h \ , the errors resulting from separate fits were greater than 
that obtained with the simultaneous fit for all values of to that we tried. 



46 



in \Vub\- For example, we can evaluate the importance of the low experimental points to 
the extraction of \Vub\ by removing them from the combined 2- fit. Including only the three 
experimental data points with > 18 GeV^, we find a consistent value of \Vub\ with only a 
~ 1% larger uncertainty. Similarly, we can evaluate the importance of having many lattice 
data points, rather than only a single point, by using only the most precise lattice point 
with a total error of ~ 9%. This allows the form factor shape to be completely determined 
by the experimental data. We find a consistent value of \Vub\ but with an even larger error 
of ~ 13%. We therefore conclude that combining all of the numerical lattice data with 
all of the experimentally- measured BABAR data minimizes the total uncertainty in \Vub\- 
Because the small error in our final determination of \Vub\ is primarily due to the power of 
the combined z-fit method, one could easily use the procedure outlined in this section to 
improve the exclusive determination of \Vub\ from existing lattice QCD calculations of the 
B iriu form factor such as that by the HPQCD Collaboration [33]. 



VI. RESULTS AND CONCLUSIONS 



Combining our latest unquenched lattice calculation of the B ntv form factor with 
the 12-bin BABAR experimental data, we find the following model-independent value for 

\Vub\:' 

X 10^ = 3.38 ±0.36. (78) 

The total error is ~ 11%, and it is nontrivial to separate the error precisely into contributions 
from statistical, systematic, and experimental uncertainty because of the combined 2-fit 
procedure used. If we assume, however, that the error in \Vub\ is dominated by the most 
precisely determined lattice point (which is not quite true, as shown in the previous section), 
we can estimate that the contributions are roughly equally divided as ~ 6% lattice statistical, 
~ 6% lattice systematic, and ~ 6% experimental. 

Our result is consistent with, although slightly lower than, our earlier, preliminary de- 

^ At three conferences during Summer 2008 wc presented a version of this model-independent analysis with 
a numerical value for \Vub\ that is l-c lower than that given here in Eq. (79). We have since improved 
several aspects of the lattice calculation, most notably reducing the statistical errors that enter the chiral 
and continuum extrapolations of f± and f\\ and, hence, /+. Equation (78) is our final result for \Vub\ 
based on the lattice data from the ensembles in Table I and the methodology of Sees. Ill and V. 
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termination of \Vub\- The reduction in central value is primarily due to a change in the 
lattice determination of the form factor, not the procedure used to determine \Vub\- Be- 
cause our new analysis uses a second lattice spacing, we are able to take the continuum 
limit of the form factor. We find that the continuum extrapolation increases the overall 
normalization of and hence decreases the value of \Vub\- Our errors are smaller than 

those of previous exclusive determinations primarily because we have reduced the size of the 
discretization errors, which are significantly smaller than in the previous Fermilab-MILC 
calculation (~ 7% — ^ 3%) because of the additional finer lattice spacing. 

Our new result is ~ l-2cr lower than most inclusive determinations of \Vub\, which typi- 
cally range from 4.0 — 4.5 x 10""^ [32]. Much of the variation among the inclusive values is 
due to the choice of input parameters — in particular that of the 6-quark mass [94]. The 
recent determination of by Kiihn, Steinhauser, and Sturm using experimental data for 
the cross section for e~^e~ — * hadrons in the bottom threshold region yields the value of mj, 
to percent-level accuracy [95], and is consistent with the PDG average [13]. Neubert has 
shown, however, that an updated extraction of from fits to B Xciiy moments using 
only the theoretically cleanest channels (excluding b — > Xg'y) results in a larger 6-quark mass 
and hence smaller inclusive value of \Vub\, thereby reducing the tension between inclusive 
and exclusive determinations [96]. 

Our result is consistent with the currently preferred values for \Vub\ determined by the 
global CKM unitarity triangle analyses of the CKMfitter Collaboration, \Vub\ x 10^ = 
S.Utoij [97], and UTFit Collaboration, |Kb| x 10^ = 3.48 ± 0.16, [98]. Further reduc- 
tion in the errors is therefore essential for a more stringent test of the CKM framework and 
a more sensitive probe of physics beyond the Standard Model. 

The dominant uncertainty in our lattice calculation of the B — > iriu form factor comes 
from the statistical errors in the 2-point and 3-point correlations. This error can be reduced 
in a straightforward manner with use of an improved source for the pion and/or additional 
gauge configurations. The statistical errors in the nonperturbative renormalization factors 
Zy and Zy can be brought to below a percent in the same way. The chiral-continuum 
extrapolation error, which is inextricably linked to the statistical errors in the correlation 
functions, can also be improved by simulating at more light quark masses and an additional 
finer lattice spacing of a ~ 0.06 fm. Presumably a better constrained chiral and continuum 
extrapolation will reduce the size of other g^-dependent errors such as those from gB'Bn, 'T'i, 
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and the light quark masses by some unknown amount as welL Use of a finer lattice with 
a ~ 0.06 fm will further decrease the momentum- dependent and heavy-quark discretization 
errors, which we now estimate with power-counting. The extraction of \Vuh\ can also be 
improved by including more experimental measurements of the B iriu branching fraction. 
This, however, will require understanding the correlations among the various systematic 
uncertainties. Given these refinements of the current calculation, an even more precise, 
model-independent value of \Vub\ can be obtained in the near future. 
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APPENDIX A: ESTIMATE OF HEAVY QUARK DISCRETIZATION ERRORS 

In this Appendix we collect the short-distance functions fi used to estimate the heavy- 
quark discretization effects. For more background, see Refs. [36-39, 83, 84]. 

1. 0{a^) errors 

We start with these because explicit expressions for the functions fi{moa) are available. 
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a. O(a^) errors from the Lagrangian 



There are two bilinears, hD ■ Eh and hiT, ■ [D x E]h, and many four-quark operators. At 
tree level the coefficients of all four-quark operators vanish and the coefficients of the two 
bilinears are the same. The mismatch function is given by 

1 1 



/ij(moa) 



2(2m2a)2' 

Using explicit expressions for l/m2 [36] and '\-/'m\ [^4], one finds 



fE{moa) 



cs(l + moa) 



moa{2 + moa){l + rrioa) 4(1 + moa)^ 
We use C£; = 1 in our numerical simulations. 



(Al) 



(A2) 



b. 0(0^) errors from the current 



There are three terms with non-zero coefficients, qTD^h, qTiT, ■ Bh, and qTa ■ Eh, 
which can be deduced from Eq. (A17) of Ref. [36]. Their coefficients can be read off from 
Eqs. (A19) [36]. When cb = Ts the ffist two share the same coefficient: 
1 Cdii)- + moa) 1 



fxirriQa) 



Im^y^a? moa(2 + moa) 2(2m2a)" 



+ 



1 



1 



(2 -|- moa)(l -|- moa) 2(1 -|- moa) 4(1 -|- moa)^ (2 -|- moa)^ 

2 



1 



(A3) 



2(1 + moa) \2(2 + moa)(l + moa) 

where the last term on the second line comes from using the tree- level di. For the third 
operator, qTa ■ Eh, 



frimoa) 



di ({1 - ce){1 + moa) 



m2a moa(2 + moa) 
-|- 4moa + (moa)^ 



4(1 + moa)2(2 + moa)2' 
where the last line reflects the choices made for ce and di. 



(A4) 



2. 0{asa) errors 

Because we improve both the action and current, the mismatch functions fi{moa) start 
at order a^, and we do not have explicit expressions for them. (The calculation of these 
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functions would be needed to match at the one-loop level.) So we shall take unimproved tree- 
level coefficients as a guide to the combinatoric factors and consider asymptotic behavior in 
the limits nioa — 0, oo. 

a. 0{asa) errors from the Lagrangian 

There are two bilinears, the kinetic energy hD'^h and the chromomagnetic moment hiT. ■ 
Bh. There is no mismatch in the coefficient of the kinetic energy, by assumption, since 
we identify the kinetic mass with the heavy-quark mass. This tuning is imperfect, but the 
associated error is budgeted in Sec. IV E. 

At the tree level the chromomagnetic mismatch is 

^''("""^ - w^y 

This has the right asymptotic behavior in both limits, so our ansatz for the one-loop mis- 
match function is simply 

ft\moa) = (A6) 
2(1 + moa) 

and errors is this function multiplied by aA. 

h. 0{asa) errors from the current 

There is only one correction at tree level, but more generally there are two for the temporal 
current and four for the spatial current. (See Eqs. (2.27)-(2.32) of Ref. [38].) 
The tree- level mismatch function ends up being the same as di. 

ff^ imoa) = -. (AT) 

■'^ ^ ^ ' 2(2 + moa)(l + moa) ^ ' 

It is, however, an accident that it vanishes as nioa 0. Therefore, we instead take 
which has the right asymptotic behavior. 
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TABLE III: Relative error from mismatches in the heavy quark Lagrangian and current for the 
bottom quark with A = 700 MeV. To obtain the totals given in the text E and X are counted 
twice, and 3 is counted twice for and four times for fj_. Entries are in per cent. 

a (fm) av{q*) niQa B 3 E X Y 
0.09 0.33 2.018 1.76 1.32 0.28 0.80 0.24 
0.12 0.41 2.617 2.48 1.94 0.39 1.26 0.33 



3. Numerical estimates 

The relative errors due to mismatches in the heavy quark Lagrangian and current on the 
MILC coarse and fine ensembles are tabulated in Table III. At the fine lattice spacing we 
take the typical av{q*) to be |, and we use one- loop running to obtain av{q*) at the coarse 
lattice spacing. The contribution errory from the a ■ E error in the current is so small 
both because c^; = 1 in our simulation and because di is small. Adding the individual errors 
given in Table III in quadrature, and taking into account multiple contributions of the same 
size, we find the total error to be 2.84% (4.16%) for /y and 3.40% (4.98%) for f± on the fine 
(coarse) lattices. 
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APPENDIX B: STATISTICAL AND SYSTEMATIC ERROR MATRICES 



In this Appendix we present the normahzed statistical and systematic bootstrap correla- 
tion matrices for the B iriu form factor, that were used in our model-independent 
determination of \Vub\- In order to facilitate the use of our result, we also show the resulting 
total covariance matrix. 
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TABLE IV: Normalized statistical (upper) and systematic (lower) bootstrap correlation matrices 
for the B iriu form factor, /_|_(g^). These should be combined with the values of f+^q"^) presented 
in Table II to reconstruct the full correlation matrices. 



g2 (GeV2 



26.5 
25.7 
25.0 
24.3 
23.5 
22.8 
22.1 
21.3 
20.6 
19.8 
19.1 
18.4 



26.5 25.7 25.0 24.3 23.5 22.8 22.1 21.3 20.6 19.8 19.1 18.4 



1.00 0.99 
0.99 1.00 
0.97 0.99 
0.88 0.92 
0.66 0.73 
0.29 0.38 
-0.01 0.07 
-0.16 -0.09 
-0.19 -0.13 
-0.14 -0.09 
-0.04 0.00 
0.05 0.08 



0.97 0.88 

0.99 0.92 

1.00 0.97 

0.97 1.00 

0.82 0.93 

0.51 0.69 

0.21 0.42 

0.04 0.25 

-0.02 0.18 

0.01 0.17 

0.07 0.20 

0.13 0.21 



0.66 0.29 -0.01 

0.73 0.38 0.07 

0.82 0.51 0.21 

0.93 0.69 0.42 

1.00 0.91 0.72 

0.91 1.00 0.94 

0.72 0.94 1.00 

0.56 0.85 0.98 

0.48 0.77 0.92 

0.43 0.69 0.84 

0.39 0.58 0.69 

0.32 0.42 0.48 



-0.16 -0.19 -0.14 
-0.09 -0.13 -0.09 
0.04 -0.02 0.01 
0.25 0.18 0.17 
0.56 0.48 0.43 
0.85 0.77 0.69 
0.98 0.92 0.84 
1.00 0.98 0.92 
0.98 1.00 0.97 
0.92 0.97 1.00 
0.78 0.86 0.95 
0.55 0.66 0.80 



-0.04 0.05 

0.00 0.08 

0.07 0.13 

0.20 0.21 

0.39 0.32 

0.58 0.42 

0.69 0.48 

0.78 0.55 

0.86 0.66 

0.95 0.80 

1.00 0.94 

0.94 1.00 



g2 (GeV^ 



26.5 25.7 25.0 24.3 23.5 22.8 22.1 21.3 20.6 19.8 19.1 18.4 



26.5 
25.7 
25.0 
24.3 
23.5 
22.8 
22.1 
21.3 
20.6 
19.8 
19.1 
18.4 



1.00 0.98 

0.98 1.00 

0.95 0.98 

0.89 0.92 

0.85 0.87 

0.88 0.87 

0.90 0.88 

0.91 0.89 

0.91 0.90 

0.92 0.91 

0.91 0.9 

0.90 0.88 



0.95 0.89 
0.98 0.92 
1.00 0.97 
0.97 1.00 
0.94 0.99 
0.94 0.99 
0.94 0.98 
0.94 0.98 
0.95 0.98 
0.95 0.98 
0.95 0.98 
0.94 0.97 



0.85 0.88 
0.87 0.87 
0.94 0.94 
0.99 0.99 
1.00 0.99 
0.99 1.00 
0.99 1.00 
0.98 1.00 
0.98 0.99 
0.98 0.99 
0.98 0.99 
0.98 0.99 



0.9 0.91 
0.88 0.89 
0.94 0.94 
0.98 0.98 
0.99 0.98 

1.0 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 



0.91 0.92 
0.90 0.91 
0.95 0.95 
0.98 0.98 
0.98 0.98 
0.99 0.99 
1.00 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 
1.00 1.00 



0.91 0.90 

0.9 0.88 

0.95 0.94 

0.98 0.97 

0.98 0.98 

0.99 0.99 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 

1.00 1.00 
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TABLE V: Total bootstrap covariance matrix for the B irliy form factor, f\{q^), derived by 
adding the statistical and systematic errors in quadrature. The elements of the matrix are given 
by Mij = o'j^(g2) X (Tj^(g2), where C/^(g2) is the total uncertainty in f+{q'^) in the z'th bin. 



g2 (GeV2) 


26.5 25.7 25.0 24.3 23.5 22.8 22.1 21.3 20.6 19.8 19.1 18.4 


26.5 


5.13 2.74 1.49 0.79 0.39 0.18 0.06 0.01 -0.0 0.01 0.03 0.05 


25.7 


2.74 1.48 0.82 0.45 0.24 0.12 0.05 0.02 0.01 0.02 0.03 0.04 


25.0 


1.49 0.82 0.47 0.27 0.15 0.09 0.05 0.03 0.02 0.02 0.02 0.03 


24.3 


0.79 0.45 0.27 0.17 0.11 0.07 0.05 0.03 0.03 0.02 0.02 0.02 


23.5 


0.39 0.24 0.15 0.11 0.08 0.06 0.04 0.04 0.03 0.03 0.02 0.02 


22.8 


0.18 0.12 0.09 0.07 0.06 0.05 0.04 0.04 0.03 0.03 0.02 0.02 


22.1 


0.06 0.05 0.05 0.05 0.04 0.04 0.04 0.03 0.03 0.03 0.02 0.02 


21.3 


0.01 0.02 0.03 0.03 0.04 0.04 0.03 0.03 0.03 0.02 0.02 0.02 


20.6 


0.00 0.01 0.02 0.03 0.03 0.03 0.03 0.03 0.03 0.02 0.02 0.02 


19.8 


0.01 0.02 0.02 0.02 0.03 0.03 0.03 0.02 0.02 0.02 0.02 0.02 


19.1 


0.03 0.03 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 


18.4 


0.05 0.04 0.03 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 
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